TRI-PP-95-53 
PSI-PR-95-14 



First Order Variational Calculation of Form Factor in a Scalar 

Nucleon— Meson Theory 



A. W. Schreiber 

TRIUMF, 4004 Wesbrook Mall, Vancouver B. C, Canada V6T 2A3 

R. Rosenfelder 

Paul Scherrer Institute, CH-5232 Villigen PSI, Switzerland 



Abstract 

We investigate a relativistic quantum field theory in the particle representation using a 
non-perturbative variational technique. The theory is that of two massive scalar particles, 
'nucleons' and 'mesons', interacting via a Yukawa coupling. We calculate the general 
Euclidean Green function involving two external nucleons and an arbitrary number of 
external mesons in the quenched approximation for the nucleons. The non-perturbative 
renormalization and truncation is done in a consistent manner and results in the same 
variational functional independent of the number of external mesons. We check that 
the calculation agrees with one- loop perturbation theory for small couplings. As an 
illustration the special case of meson absorption on the nucleon is considered in detail. 
We derive the radius of the dressed particle and numerically investigate the vertex function 
after analytic continuation to Minkowski space. 



PACS numbers: ll.80.Fv, 24.10.Jv 
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1 Introduction 



There are many instances in relativistic quantum field theory where perturbation theory is 
inadequate, in particular when the coupling constant characterizing the interaction becomes 
too large. An alternative approach, which remains valid in the strong coupling regime, is the 
use of a variational principle. Although there exists a sizable body of literature on the use of 
variational methods in relativistic quantum field theories |^, ^, ^, |^ it is fair to say that 
these methods are less widely used in this area of physics than the standard Rayleigh-Ritz 
variational method in non-relativistic quantum mechanics. This is largely due to the fact that 
one has to deal with a system with infinitely many degrees of freedom, the requirements of 
relativistic invariance and renormalizability of the theory. 

The problem of the infinitely many degrees of freedom, however, is also present in non- 
relativistic field theories and here there exists a variational method which deals with this 
problem very successfully - Feynman's variational treatment of the polaron |^, |^, [1^, 11 



A very early attempt to use this method in a relativistic setting was made in Ref. []T3| and 
was recently reinvented and extended by us fl^, [1^, |16|. (The method has also previously been 
adapted for use in non-relativistic nuclear physics [l^, |18|, |l^). The particular Lagrangian which 
was considered in these papers was that of two scalar particles, one heavy (henceforth referred 
to as the 'nucleon') and one light (the 'meson'), interacting via a Yukawa coupling. This variant 
of the Wick-Cutkosky model ||20|, ^ was chosen because of its simplicity, its similarity to several 
realistic Lagrangians (such as those of QED and nucleon-meson theories) and because it has 
traditionally been used for the study of non-perturbative techniques [^, |2^, ^ |2B|, 

The quantities of interest in a relativistic field theory are the Green functions of the theory. 
The work in Refs. |T^, [T^ concerned itself with the calculation of a subset of these Green 
functions in the quenched approximation for the nucleon. In particular. Green functions with 
2 external (on-shell) nucleon legs and any number of external meson legs were calculated in 
Ref. [ig. 

Feynman's variational method involves the use of the path integral formalism. As is well 
known, in general it is only possible to do gaussian path integrals analytically. Because of 
this he made use of a particular approximation - Jensen's inequality which can be exploited 
for variational purposes. Equivalently, this corresponds to the first order term in a cumulant 
expansion (see, for example, Ref. |29|). In Refs. |1^, [1^, the variational calculation of the 



nucleon propagator was performed using the same approximation. The only quantity of interest 
there was the residue of this propagator, the position of the pole being defined as the (non- 
perturbatively) renormalized mass of the nucleon. (In addition, the 'width' of the nucleon 
was calculated - this quantity is not zero as the Wick-Cutkosky model is inherently unstable.) 



In Ref. 0] the general (2 + n)-point function was then written down to zeroth order in the 



cumulant expansion. Although already qualitatively interesting at this order (for example, even 
at zeroth order these Green functions contain parts of all Feynman diagrams to all orders in the 
coupling), these Green functions clearly had some unwanted numerica/ features. In particular, 
at this order all dependence on the external momenta is gaussian, resulting in, for example, 
a gaussian form factor. Furthermore, one of the strong points of Feynman's treatment of the 
polaron was that agreement with perturbation theory was assured to the order in the cumulant 



expansion in which one calculates. Hence, the (2 + n)-point functions calculated in Ref. |jT6 
only agree with perturbation theory at tree level. 
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In the present paper we calculate the (2 + n)-point function to first order in the cumulant 
expansion, therefore greatly improving the numerical results. More importantly, it will be 
shown that the variational equations for any of these Green functions are in fact the same, 
i.e. once one has determined the numerical value of the variational parameters from one of 
these functions, say the propagator of the nucleon (n = 0), the other ones may be calculated 
with the use of the same parameters. It will also be shown that it is possible to truncate the 
(2 + n)-point function consistently; that is, the propagators on the external nucleon legs have 
poles at the physical nucleon mass, as they should. 

An important ingredient of the method has been the use of the particle representation 
of a field theory ^ |3^, which has undergone something of a revival in recent 



years p2|, |3^, ^ |3^, We start in Section || by briefly outlining how one derives 

the expression for the quenched (2 + n)-point function in the particle representation. Having 
done this, we perform the actual variational calculation in Section ^, resulting in the 'master 
formula' for the (2 + n)-point function at first order in the cumulant expansion. An important 
check, that this formula agrees with perturbation theory to one-loop order, is relegated to the 
Appendix. It should be noted that throughout this derivation we work in Euclidean space- 
time. The Green functions can however be continued to Minkowski space. This is done, for the 
(2 + l)-point function, in Section ^ and numerical results are derived and discussed in Section 
We conclude with some general observations in Section H 



2 The (2+n)-point function in the particle representation 

Before we turn to the calculation of the general quenched Green function with two external 
nucleon and n external meson legs in first order in the cumulant expansion, we need to derive the 
particle representation of this function. For the propagator this was already done in Ref. |]14| , 
so here we shall limit ourselves to briefly summarizing the relevant steps involved. We shall 



follow the notation of Ref. ||Tj], and the reader who is interested in the details is referred to 



this reference. The Wick-Cutkosky Lagrangian is given, in Euclidean space-time, by 

C = \ {d,^f + Im^o^' + \ {d,^f + ^mV^ - (1) 

where Mq is the bare mass of the heavy nucleon, (p and $ are the fields of the meson and nucleon 
respectively and m is the meson mass. The latter does not get renormalized in the quenched 
approximation while the nucleon mass does. The corresponding generating functional for the 
general Green function is written in terms of the sources J and j of the external nucleons and 
mesons as 

Z [J,j] = jv^Vip exp if] + (J, $) + (j, if) ) . (2) 

Here the action is given by 

S[^,^]= j d^x C{^{x),^{x)) (3) 

and we have defined 

(J, j d^x J{x)^{x) etc. (4) 

One implements the quenched approximation by integrating out the nucleon field and setting 
the resulting determinant equal to unity. After functionally differentiating twice with respect to 
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the nucleoli source J{x) and then setting this source to zero we obtain the generating functional 
for the Green functions with two external nucleon legs and an arbitrary number of meson legs. 
It is given by 



6J{x) 6J{0) 
j Vip <x\ — 



(5) 



j=o 



> exp 



--(^,(-n + m2)y,) + (j,(^) 



This generating functional is written in terms of the meson field y?. In order to implement the 
particle representation of the theory in terms of the path of the nucleon we exponentiate the 
propagator and write the resulting relativistic matrix element in terms of a proper time path 
integral [PU|: 



< X I exp 



>= 



x(l3)=x 
a;(0)=0 



T>x{t) exp 



dr 



■ (6) 



We have given the particle the "mass" 1 which corresponds to the 'proper time gauge' Our 
variational results can be shown to be gauge-independent although some variational parameters 
are not One may now perform the functional integration over the meson fields to obtain 
the generating functional of connected Green functions 



^conn b> = const f dp CXp [-^^0 



x{l3)=x 
x(0)=0 



Vx{t) exp(-S'cfr [x(r), j] 



(7) 



Here the effective action is given by 

'S'eff [xir)J] = S [x{t)] - g J d yj{y) dr < y\ — ^ 



x{t) > 



with 



S[xiT)] = So[xiT)] + Si[xir)] 

rP 1 

50 [x{t)] = dr -x^ 

Jo 2 

g^ fl3 rP 

51 \x{t)] = / dti / dt2 < x{ti) 

2 Jo Jo 



X{t2) > 



(8) 

(9) 

(10) 



Finally, one may obtain the Green function with n external meson legs by differentiating n times 
with respect to the source j and then setting this source to zero. After Fourier transforming 
to momentum space and removing the external meson legs as well as the overall momentum 
conserving delta function one obtains 



G2nip,p'', {q}) = const. / d(3exp — ^Mq / d'^x e *^ ' 

JO 2 J 



x{f3)=x " 

Vx{r) n 

x{0)=0 -^-^ 



-S[x{r)] 



(12) 
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Here the definition of tlie meson momenta {q} and of the final nucleon momentum p' is such 
that the latter is in-going while the former are out-going. The integration parameters Xj have 
an obvious physical significance - they mark the (proper) times at which the external mesons 
interact with the (bare) nucleon. 

Note that the external nucleon legs are still contained in this expression. As the nucleon 
mass needs to be renormalized their truncation is not as trivial as that of the meson legs. In 
particular, the nucleon mass is a dynamical quantity which needs to be calculated in terms of the 
bare nucleon mass. Of course, if one could calculate the Green functions of the theory exactly, 
the poles of the external nucleon legs would all be automatically at the same position as the pole 
of the propagator. In an approximate treatment, however, this is not automatically fulfilled. 
In fact, it is non-trivial to even isolate these poles, let alone to have them occur at the correct 
position. As will be seen in the next Section, it is indeed possible to perform the truncation 
of the general Green function in a consistent manner within the context of the variational 
calculation described here. The result (see (Eq. ^) ) will be a rather compact expression for 
the general truncated (2 + n) - point function, which encapsulates, as we shall see, information 
contained in any order of perturbation theory as well as genuinely non-perturbative effects. 

3 Variational Calculation of the (2 + n)-point function to 
first order in the cumulant expansion 

The variational treatment is based on the decomposition of the action into a trial action St 
plus a remainder 

S = St + S - St = St + AS (13) 
and uses the cumulant expansion to first order (which is related to Jensen's inequality): 

< e~^^ ^ e-<^^>^* . (14) 

Here the action S[x] is related to S[x] through 

n 

S[x] = S[x{t)] +ip' ■ X - i^Qi- x{Ti) (15) 



=1 



and thus contains all the dependence on the nucleon's path in Eq. (|T2|). The first order cumulant 
has the particular property that it is stationary (for real actions even minimal) under arbitrary 
variations of 5*^. We shall perform the variational calculation in momentum space so that the 
averages [| in Eq. (|14D are defined by |14 



JViO exp(-St) , , 

< C >^ = - ' . 16 



^ These averages in momentum space were denoted by << .. >> in Rcf. p4[ | in order to distinguish them 
from averages in coordinate space which were indicated by a single bracket < .. >. In the present paper we 
shall only deal with momentum space averaging, so as no confusion can arise we shall simply make use of the 
single brackets. 
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a 
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Figure 1: The definition of tlie relative times and momenta. Note tliat tq = 0, r„_|_i = /5, 
Pi = -p and pn+i = p'. 



wliere 



rx(f3)=x 

d X T>x 

Jx{0)=0 



Vi = I d^x I Vx (17) 
To tliis order in tfie cumulant expansion tlie (2 + r;,)-point function is given by 



G2,n(p,p'; {q}) = const. / d(3 

Jo 



U9 

i=l 



/3 



exp 



:i8) 



Tlie actions S and St are expressed in terms of the absolute times Xj. It will at times turn out 
to be more useful to define relative times a^, as indicated in Fig. [l|. It should be noted that 
Fig. |1| is not to be understood as a Feynman diagram - it only serves to label the interaction 
times of the external mesons with the nucleon. Furthermore, in this Figure we have ordered 
these interaction times such that n < r2.... < r„. It is always possible to do this, provided we 
sum over all permutations of the g^'s afterwards. 

We shall do the path integrals which are involved by expanding the general path x{t) into 
its Fourier components: 



x[T) = X — + ^ — — Ok sm 



k=l 



The action S now becomes 

S = S[x{T)]+t 



kn 



n+l 



i=l 



k=l 



(19) 



(20) 
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where 



in) 



2v^A . fknr- 



(21) 



X 



'2/3 60 (22) 

and the definition of the momenta Pi is indicated in Fig. |I[ 

The question now arises as to what the form of the trial action St should be. If one were to 
calculate all the corrections to Eq. (0) the functional form of this trial action would in principle 
be arbitrary, as would be the actual numerical values of the variational parameters. In practice, 
of course, one does not calculate the corrections, except the lowest order one which is needed 
to make the variational principle work. It is this fact that makes the form of the variational 
action assume importance. Indeed, for best results, the trial action should encapsulate as much 
of the information contained in the true action as possible. In particular, in order to reproduce 
the lowest non-trivial order of perturbation theory, the trial action should reduce to the free 
action in the small coupling limit [Q. On the other hand, it should be of such a form that 
the mathematical manipulations required become possible without further approximations, as 
otherwise one looses control over the variational properties. In the present case this simply 
means that the trial action must be quadratic in the Fourier coefficients bk so that the path 
integrals may be performed analytically. One is therefore restricted to the following general 
form: 



St = St + i 



n+l 



'-,po bo 



J2 ^iPi -J2pk Qk ■ h 
i=l k=l 



(23) 



where St is given by (as in Ref. [T^ 



St 



00 



k=0 



(24) 



Note that in the limit in which the coupling goes to zero the true action S becomes 



k=0 



1 2 n+l 

-60 ■ ^ QiPi 

y i=i 



00 

E 

k=l 



(25) 



so that Ak and pk must go to 1 in this limit. 

The above trial action ensures that upon expansion in the coupling constant, and without 
having to solve the variational equations, the one-loop perturbative result is obtained. However, 
as indicated in Section ^, we also demand that for arbitrary coupling one should be able to 
truncate the general Green function consistently. It is not obvious that this is possible in 
general. Indeed, we have not been able to perform this truncation with the simplest choice of 
pk being a constant independent of k. On the other hand, if one chooses pk to be proportional 
to the 'profile function' Ak 

Pk = \Ak , (26) 

a consistent truncation becomes possible. Here A is a variational parameter which will turn 
out to be precisely the same quantity which already appeared in the study of the propagator 
in Refs. E 



^. Actually for n = Eq. (p3|) with the choice ( |26| ) is just the ansatz which 
was used there for 'momentum averaging' of the 2-point function. However, for this special 
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case it is more a convenience than a necessity to write the variational parameter po as XAq 
whereas for n > the /c-dependence chosen in Eq. (p6D is essential. At the present stage we 
have not investigated other parameterizations of the trial action in a systematic manner. As 
the parameterization of the trial action is connected to the analytic structure of the resulting 
Green functions it would clearly be a rather interesting exercise to determine what the most 
general allowable is. 

One may now calculate the various quantities required for the evaluation of the Green 
functions in Eq. (p!8[). Writing 

S = So + Si (27) 
one obtains, up to irrelevant overall constants, 



/ 



Vxe 



-St 



■ oo 1 ■ 
.k=0 ^k. 



^ fc=i 



(2^ 



and 



< Sq — St >St~ ~7r 



p2 



A-y(l + Ao) 



°° / 1 



k=0 



oo g(n) 



E 

k=l 



2 ^ 



A-y(l + ^.) 



(29) 



For brevity we have defined the quantity P to be 



n+l 
i=l 



(30) 



Combining these terms yields 



exp 



< 5o - > 



Vx e 



exp 



oo . ^ 



A(A-2)"+l 2 
^ E ^iPi 



i=l 



(31) 



The remaining quantity which needs to be calculated is < Si > g^. It is given by 
^2 r d\' 1 



< ^1 >st= -- 



9 



2 Jo 



dtidt2 



(27r)4 g'2 ^ ^2 



< exp[zg'-(x(ti)-x(t2))] >5, • (32) 



The path integrals in this expression can be performed analogously to the equivalent expression 
for the propagator in Ref. |T3. Furthermore, the integral over q' may also be done with the 
help of the momentum space representation for the propagator 



1 



g/2 _|_ 



_ 1 

2 ~ 2 



du exp 



u 



(33) 



One obtains, after shifting the u integration in the same way as in Ref. 

.2 in' 



< Si >si- 



9 



dtidt2 — 



/i2(a,T) Jo 



du e \ mfi{(T,T), , m 1 . (34) 
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The quantity jl {a,T) is defined in terms of tlie pseudotime /j, {a,T) as in Ref. fl^, i.e. 

+ /3E>- (35) 



Here 

A 

Afc — 1 

KTT 



Ao(3 ^ Ak 



( k7Tti\ ( kTTt2' 



(36) 



and we liave introduced tlie relative and total time a = ti — t2 and T = (ti + 12)/2. Finally, 
the function e{s,t,u) is defined as 

e(s, t, m) = exp I - — — Y " j • (37) 



It should be noted that the interaction term Eq. ( p4D looks identical to the interaction term 
written down for the propagator in Ref. []T3[, Eq. (97). Indeed, all the dependence on the 
number of external mesons is hidden in the (four dimensional) quantity W^'^^: 



iy(«) = ^p-^^EA.Ql"^ (38) 

For n = this just reduces to a times the nucleon momentum p, while for general n it contains 
the full complexity necessary for the (2 + n)-point function G2,n- Indeed, using the standard 
formula for the summation of a series of cosines weighted by l/k"^ [^, we may evaluate the 
sums appearing in W^'^^ to obtain 

in 

W^^"^ = + Pn+l) + 1^ E 9^ (1^^ - ^1 1 - - ) . (39) 

^ 1=1 

Hence W^""^ depends on the ordering of the proper times Ti at which the nucleon interacts with 
the external mesons relative to the proper times ti^2 at which the internal meson is emitted 
and absorbed. The form of Eqs. (^) and (|39|) bears some similarity to the "phase averaged" 
approximation for the Green function and corrections to it which have been investigated in 



Ref. 42 



It is instructive to eliminate the absolute values appearing in Eq. (p9| ) by explicitly breaking 
up the sum into three contributions: one contribution (z = 1 to b, say) where the internal 
meson's interaction times ti^2 are both larger than all t^, one contribution where one of them 
is smaller and the other one is larger than all (z = 6 + 1 to a — 1) and one contribution 
where they are both smaller than all {i = a to n). This division of the sum is illustrated in 
Fig. 1^. Similarly to Fig. |I], Fig. ^ should not be interpreted in the sense of a Feynman diagram, 
although it is closely related. It is a diagrammatic representation of the (in this case averaged) 
action (not an amplitude) in the particle representation. It serves to indicate the relative 
ordering of the proper times involved. Furthermore, it also defines the momenta: similarly to 
a Feynman diagram, four-momentum is conserved at each vertex, with the proviso that one 
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b+1 





a+1 ♦ 



b+1 



a-1 



Figure 2: The relative ordering of the interaction times of the external mesons with respect to 
those of the internal one. It is assumed that a > 0. 



associates four-momentum with the internal meson line (the momentum flowing though the 
loop has already been integrated over). 
One obtains 



= sign(cr)i [rfe+i - min(ti,t2)K+i + '^iPi 

^ i=b+2 

+ [max(ti,t2) - ra-i]pa | a> b + 1 



O-Pa 



a = b+l. (40) 



(For the case a = b + 2, the sum appearing in Eq. (^) is defined to be empty.) In other words, 
W^"^^ is equal to (up to a sign, which is not relevant as only the square of W^'^^ enters) the 
integral of the proper time multiplied by the nucleon's four-momentum for the duration of the 
exchange of the internal meson. The reason for the appearance of this quantity will become 
clear once we make the connection to perturbation theory. Particular cases of interest, for the 
purpose of the eventual truncation of the Green function, are those where ti^2 are both smaller 
or both larger than all Tj, in which case W^"'^ = api and apn+i, respectively. 

As was shown in Ref. [jl4|, if one is interested in the Green function for the case where 
the external nucleon legs {not the internal ones) are on-shell, then one should consider the 
above expression in the limit where the proper times corresponding to the external legs (i.e. ai 
and On+i, and hence /?) tend towards infinity. To implement this limit, the following equality 
derivable from the Poisson summation formula [^, ^ is extremely useful: If F{k) is only a 
function of kir/jS and if it is even, then up to exponentially small terms in (3 



y F{k) ~ ^ / dEF 



(41) 
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As was shown in Ref.|jT^, Ak satisfies these criteria, which results in the following simplifica- 
tions: 

hm /i^ a, T ^^\o) = - \ dE — — "^—^ . 42 

P^oo TT Jo A{L) E^ 



and 



oo 

2 lim y 



\ogAk + 



A, 



21 

IT Jo 



dE 



log A{E) + 



A{E) 



- 1 



logA(O) + 



+ logA(O) 
1 



1 



1 



A{0) 



AiO) 

n+l 
i=l 



(43) 



Here Q is identical to the kinetic term defined in Ref. 14 



Using the above simplifications for the case of on-shell external nucleon legs, as well as the 
fact that the integrand appearing in the interaction term < Si is even in a, one can now 
write down the expression for the (untruncated) (2 + n)-point function: 



G'2,n(p,p'; {q}) 



No 
29 



'n+l 



n<l^} 

■ exp 



5: n W exp {m^ + 2n+pni - (1 - A)^ 



(44) 



Here No is defined to be 



87r2 



dti / dt2 



du e ■mjj,{a), 



/x(cr) 



u 



AiO) 



exp 1 



(45) 



and it should be noted that G2,n{p,p''i {l}) is normalized in such a way that the correct tree 
level result is obtained in the limit in which the coupling vanishes: 



1 n+l 

th: E n 



2^ 



^9 



V{{q}) i-- 



1 P^'' + Mp\y, 



(46) 



Furthermore, for non-zero couplings Eq. 
come meaningful. This is because yU^(o") 
diverges logarithmically. Indeed, the mass renormalization was already performed in Ref. |jl4 



needs to be renormalized in order for it to be- 
er for small a and hence the integral in Eq. 



by demanding that the pole of the nucleon propagator occurs at the physical nucleon mass. 
Actually, the nucleon propagator is just the special case of the Green function in Eq. (0) with 
n = 0. It is therefore not surprising, but nevertheless crucial for the internal consistency of the 
variational approach described here, that with the same renormalization on the poles of the 
external legs we arrive at the same relationship between the unrenormalized nucleon mass Mq 
and the physical nucleon mass M as one did in Ref. ||14|| : 



XaM 



u 



(47) 



11 



■A 




It should be understood, of course, that in principle one should regularize both Eqs. 
and (^71). Only after substitution of the latter into the former does everything become finite 
and therefore meaningful with the regulator removed. 

Because the meaning of W^"^^ depends on the relative ordering of the proper times ti_2 with 
respect to the interaction times r,, the ti - t2 integration region in Eq. ( PI ) needs to be broken 
up in order to do explicit calculations. In particular, in order to isolate those parts of the 
expression which contribute to the dressing of the external legs one needs to separate out the 
parts of the integration region for which ti 2 are either both smaller than all Tj (Region {A) in 
Fig. 1^) or both larger than all (Region (B) in Fig. |^). Explicitly, in the limit in which ai 
and become very large, the relevant integrals for these two regions become 



p poo 

/ dtidt2 — > / daiai — a) 

J{A) Jo 

r poo 

/ dtidt2 / da(an+i — o") , (4^ 

J(B) Jo 



l{B) 

respectively. Here we have made use of the fact that in regions (A) and (B) the integrand does 
not depend on the overall time T. 

It is important to note that there is no further dependence on ai and dependence 
coming from the integration region (C). This allows one to perform the integrals over these 



two proper times. Comparing with the results for the propagator in Ref. |jT^, one sees that 
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the result of these two integrations just yield the external nucleon propagators, including their 
residues: 



^ dai expj -y 



+ 



[M'+p^] (1 
da n 



11 

Atx"^ Jo /i2((T) Jo 

3 



du I e (m/i((T), — ^-r , u 



e mjj,[a), — , u 



exp 



da 



a 



li\a) Jo 



ctM e ■mjj,[a), — — — , u 



^ 2Z/No 

p2 + M2 

Here i is 1 or n + 1, respectively, and Z is the same expression for the nucleon propagator's 
residue as derived in Ref. [|I^. 

By removing the external legs one may now arrive at the main result of this paper, the 
truncated (2 + n)-point function G^2niPyP''i valid for on-shell external nucleon momenta: 



Gnp,p';{q}) 



dai exp 



a. 




^{{M'+pfKi-ii-xr) 



(C) /i2(a) 



^— ^^a^e m/i(a), — — , u 
e m/i((T), , u 



(49) 



In order to appreciate the information contained in Eq. ( ^91) it is very instructive to expand the 
exponential in the second and third lines in something akin to a perturbative expansion (i.e. 
expand in g'^, disregarding the dependence of iVo, ^ and \ on g). The leading term contains 
the tree diagram and for this one the integrations over the proper times at may actually be 
carried out, resulting in the appropriate propagators for the internal nucleon legs. The next 
term contains diagrams occuring at one loop order. Because the proper times ti and ^2 at 
which the internal meson couples to the internal nucleon line are integrated over, they occur 
in all possible permutations with respect to the external times r^. These permutations give 
rise to all the individual Feynman diagrams at this order. Similar reasoning applies to the 
higher order terms - at order g'^ four integrations over the (two) internal mesons' interaction 
times must be carried out, each ordering with respect to the times Ti corresponding to a piece 
contained in a particular Feynman diagram, etc. Of course, there is no reason for a particular 
one of these terms to reproduce a given Feynman diagram exactly. This is only guaranteed for 
the one-loop result and is explicitly checked in the Appendix. 

Using this procedure, it becomes possible to address the question whether or not Eq. (^9]) 
explicitly contains all possible Feynman diagrams of any order . (In addition, it contains 
additional information implicitly because, after all, Nq, fi and A do depend on the coupling. A 
diagram by diagram identification of this implicit dependence on the coupling constant does 
not seem possible.) In fact, it is easy to convince oneself that Eq. (^Up does contain, in an 
approximate way, almost every diagram occurring at any order in perturbation theory, except 
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for an important exception. Because of the restriction of the ti — t2 integration to Region (C) 
in Fig. 1^, diagrams of fourth order or higher in which ti and t2 are both smaller than Ti or 
both larger than r„ (excluding, of course, the one-particle reducible diagrams taken away in 
the truncation) are not contained explicitly. It would be very interesting to see to what extent 
this feature depends on the form of the variational function pk in Eq. (p3|). 

Although Eq. ( p9D is a very compact expression for all the Green functions in the one 
nucleon sector of the theory, it is in general not possible to continue without specifying the 
number of external meson legs n. The simplest case (apart from the nucleon propagator, which 
has already been worked out in Refs. [|T^, ^]) is the vertex function, i.e. n = 1. We shall 
discuss this function in the next Section. The rather more complicated case of the (2 + 2)-point 
function ('meson - nucleon', or 'Compton' scattering) will be discussed in a future paper. 



4 The vertex function 

Setting n = 1 in the expression for the general truncated Green function G2^j^{p,p'; {q}) in 
Eq. (p9D yields the vertex function. For this function, the only possible values for a and b, as 
defined in Fig. ^, are 2 and 0, respectively. Hence W^"'^ takes on the form, for positive cr, 

= |(Pi+P2) + (T-ai)gi . (50) 

It is important to note that although W^^^ depends on ai, the integral in the last line in 
Eq. (^) does not. The ai dependence disappears by a shift of variable: the integration region 
(C) for the vertex function is given by, in the limit of ai ^ oo appropriate for on-shell external 
nucleons, 

p poo pax^u j1 

\ dtidt2 = da dT. (51) 

J{C) Jo Jai-a/2 

Letting T — > T — ai eliminates both the ai dependence in the integral limits as well as in 
W^^^ . This is absolutely necessary - after all, the variable ai has already been integrated over, 
resulting in one of the external propagators of the untruncated Green function. With this shift, 
the scattering amplitude is given by 

^2,i(p,p';g) = t du r ^J'^^'dT (52) 

Nq [StT^ jo jo /i (o-) J-cr/2 



■- m /X a ^ — u 

2 V u /i2(a) , 



■ exp 

Note the extra factor of Z relating the Green function in Eq. ( ^9]) and the corresponding 
amplitude. 

As it stands, the integrand in Eq. (^) looks as if it may diverge for small a due to the factor 
l//i^(cT). It is easy to see that this is actually not the case by scaling the T integral by a. That 
we do not need an additional (coupling constant) renormalization is, of course, a special feature 
of the Wick-Cutkosky model in the quenched approximation which is a superrenormalizable 
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theory. By also making use of the fact that pf = P2 = — M^, one obtains the form factor as a 
function of the momentum transfer q"^ = —qf in Minkowski space 



■ exp 




--\m fi (T, 
2 \ u 



Note that, in contrast to perturbation theory, for example, the change from Euclidean space to 
Minkowski space was trivial here as there are no integrations over momenta to be performed. 
Of course, Eq. (|53| ) only defines the amplitude in the region in which the integrals converge. 
This is the case for spacelike momentum transfers, as well as (unphysical) timelike momentum 
transfers as long as < AM"^. In these regions the amplitude is purely real. For timelike 
momentum transfers larger than AM"^, however, one needs to analytically continue Eq. (p3D 
before obtaining numerical results. It is possible to do this if the analytic structure of /U^(ct) 
is known. In particular, if /U^(cr) has no cuts or poles for Re(cr) > 0, then one may deform the 
contour of the a integration in Eq. (|53|) to the imaginary axis. If one gives a small positive 
imaginary part, one must perform the deformation in the positive half-plane so that one does 
not get a contribution from the contour at infinity. Setting a = is, one obtains 

A'o [ Stt^ jo jo fi^{is) 
■ exp 




It is clear that for > 4M^ the amplitude is no longer purely real. Indeed, it develops a 
cut, corresponding to the fact that this is the physical region for nucleon pair-production. 
It should of course be remembered that the calculation has been performed in the quenched 
approximation; hence the numerical results in this region certainly do not correspond to what 
one would obtain in the full theory. 



5 Numerical Results and Discussion 

Our numerical results for the vertex function will be presented as a function of the dimensionless 
coupling constant 

We only consider a below the critical coupling ac — 0.82. For couplings larger than ac the 
variational equations do not have real solutions due to the well-known instability of the Wick- 



Cutkosky model|14, 15 



We first investigate the vertex function for low values of the momentum transfer. Writing 
Eq. (ID as 

A2,i(g^) = A2,i(0)F(g2) (56) 

one may easily read off the effective coupling at zero momentum transfer and the elastic form 
factor for scattering of mesons from the dressed particle. The effective coupling takes on a 
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particularly simple form as it may be simplified through the use of the variational equation for 
the parameter A (see Ref. [jl5|, Sec. 5.1). One obtains 



2(7eff = ^,l(0) 



2l 
A 



(57) 



Like the famous anomalous magnetic moment of the electron in QED the effective coupling is 
enhanced because the 'velocity parameter' A is always less than unity. It is interesting to note 
that in zeroth variational order we had obtained the enhancement = g Aq |]16| which can be 
shown to be less than the one obtained in first order. Indeed, solving the variational equations 
perturbatively for m = we obtain 



gcs - g 
9 



2-K 



+ ... 
+ ... 



zeroth order 
first order ) . 



(5^ 



In addition Aq varies much more among the different parametrizations in Ref. [|15] of the re- 
tardation function than does the parameter A which is essential for determining the critical 
coupling. 

The elastic form factor is given by 



exp< 



9^ 



du 



da 



a 



^2(a) 



dt 







1 — exp 



exp 



2^2 



AV 



2 2/ \ -'- 

--m /i (a) 

2 u 



'l-e)q^u 



u 



(59) 



From its low-g expansion we deduce the mean square radius of the dressed particle to be 

-.2 /.l /.oo ^3 

/o 



r 



167r2 



A^ 



du u 



OO (J 

da 







^^\a) 



exp 



\ ( mV(^) — 
2 V u 



u X^M^a^ 
- + — ^—u 



/i2(a) 



(60) 



In the limit A l,/i^(a) 
result 



a we can perform the cx-integration and recover the perturbative 



pert 



47r2 



du 



u 



ujm^ 



(61) 



For large spacelike momenta <^ — the elastic form factor (|59| ) approaches the constant 
value 



exp< 



9' 



^vr^ Jo 



du 



OO (J 

da 



/i2(a) 



exp 



- m /i a + u 



(62) 



which is identical with the constant A^^i = XZ/Nq encountered in the calculation of the residue 
[]T3[ . Although the large g-limit certainly is outside the region of applicability of the quenched 
model the result still is meaningful: for asymptotically large momentum transfers the meson 
"sees" the bare particle inside the Wick-Cutkosky polaron. 

The effective coupling constant and the radius of the dressed particle are plotted as a func- 
tion of the dimensionless coupling in Figs. ^ and ^ We have employed several parametrizations 
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'Feynman' 


'improved' 


'pxtpndpd' 


'variational' 


perturbative 


0.1 


0.0480 


0.0481 


0.0481 


0.0481 


0.0470 


0.2 


0.0697 


0.0697 


0.0697 


0.0697 


0.0664 


0.3 


0.0878 


0.0879 


0.0879 


0.0879 


0.0813 


0.4 


0.1047 


0.1049 


0.1049 


0.1049 


0.0939 


0.5 


0.1216 


0.1220 


0.1220 


0.1220 


0.1050 


0.6 


0.1397 


0.1404 


0.1404 


0.1404 


0.1150 


0.7 


0.1610 


0.1622 


0.1622 


0.1623 


0.1242 


0.8 


0.1937 


0.1980 


0.1979 


0.1985 


0.1328 



Table 1: First-order root-mean-square radius (in fm) of the dressed particle from Eq. (|60|) for 
various parametrizations of the retardation function. The heading 'Feynman' gives the result in 
the Feynman parametrization whereas 'improved' refers to the 'improved' parametrization. In 
both cases the parameters have been determined from the variational calculation for the nucleon 
self-energy in Ref. ([|15|) . The 'extended' parametrization was introduced in Ref. ([|16|) for 
additional use in meson-nucleon scattering. The radius calculated with the solution of the 
variational equations is denoted by 'variational'. For comparison the perturbative result is also 
given. 



of the retardation function in the quadratic trial action 0. The nucleon mass is taken to be 
939 MeV, while the meson mass is fixed at 140 MeV. For comparison, the perturbative and 
zeroth order results are also shown. Table 1 demonstrates quantitatively that the results of 
the variational calculation are practically independent of which ansatz for the trial action is 
chosen. Indeed, if one leaves the form of the profile function free and lets the variational prin- 
ciple determine it (corresponding to the 'variational' trial action in Ref. |T^), the results for 
the effective coupling and radius are nearly indistinguishable from those of the 'improved' or 
the 'extended' parametrization. Furthermore, the results agree with perturbation theory for 
small couplings while they deviate quite strongly for large couplings. 

Even though there are large numerical differences between the perturbative and variational 
calculation of the effective coupling and radius, one does not observe a qualitative difference. 
The situation is dramatically different for the form factor away from zero momentum transfer. 
The real part is plotted, for three different couplings, as a function of the meson's momentum 
transfer squared in Fig. ^ For < 4M^ the representation of Eq. (^) was used, while 
for > 4M^ we used the analytically continued representation of Eq. (|S3). The physical 
region for meson absorption is < 0, while > 4M^ corresponds to pair-production. (The 
unphysical region < g^ < 4M^ is also shown because the matching of the two representations 
at g^ = 4M^ provides a useful check on the accuracy of the numerical integrations.) 

In general, the perturbative and variational calculations tend to approach each other at 

^See Refs. |l^ for more details about these trial actions; the 'Feynman' trial action is of the same form 
as that used in the variational calculation for the polaron |6| while the 'improved' trial action has the same 
ultraviolet behavior as the true action. The 'extended' parametrization gives in addition the correct threshold 
behaviour in zeroth order meson-nucleon scattering and leads to the lowest value of the variational functional 
among the three parametrizations. 
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1.6 




1.6 




Figure 4: The effective coupiing strengtli as a function of tlie dimensionless coupling a : (a) in 
zerotfi order, (b) in first order of tfie variational calculation. The solid (dashed) line corresponds 
to the use of the 'improved' ('Feynman') trial actions, while the dash-dotted line is the result 
from the 'extended' parametrization. The dotted line denotes the perturbative result. 
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0.20 



I — I 




0.0 0.2 0.4 0.6 0.8 1.0 

a 

Figure 5: The mean square radius as a function of the dimensionless couphng a : (a) in 
zeroth order, (b) in first order of the variational calculation. The various curves have the same 
meaning as those in Fig. ^. 
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Figure 6: The real part of the amphtude as a function of the momentum transfer squared, for 
from top to bottom, a = 0.2 (a), 0.4 (b) and 0.8 (c) The sohd hne corresponds to the use of 
the 'improved' trial action, while the dotted line is the perturbative result 



20 



large values of |g^|. In the area just above threshold, however, major differences arise as the 
coupling is increased. Whereas the perturbative amplitude increases smoothly in magnitude 
as the coupling is increased from a = 0.2 to a = 0.8, the variational amplitude decreases 
from a = 0.2 to a = 0.4 and actually becomes negative for a = 0.8. As the amplitude has 
to be positive at threshold ( Eq. (^3])) is always positive) this necessitates a rapid oscillation 
very close to threshold. Furthermore, near threshold the variational amplitude is numerically 
greatly enhanced (at a = 0.8 by orders of magnitude!) as compared to the perturbative result. 

A somewhat similar behavior is exhibited by the imaginary part of the amplitude, shown in 
Fig. 1^. This is, of course, only non-zero above threshold. Again the variational calculation and 
perturbation theory approach each other for large and the two calculations are qualitatively 
similar for the relatively small coupling of a = 0.2 ( Numerically, though, there is already 
a factor of two between the two calculations at the peak.) In contrast to the real part of 
the amplitude, the imaginary part of the variational amplitude is still increasing at modest 
couplings (a = 0.4), but again at a = 0.8 large oscillations show up. 

The pronounced quantitative and qualitative differences between the leading order pertur- 
bative calculation and the variational one are due to the fact that the latter has not been 
expanded in powers of the coupling constant. Let us writes Eq. (p3) as 



^2i(?^) = const, exp I — Re H(g^, a) + z Im E!(g^, a) i , (63) 

1 27r L -I J 

where the dependence of ^(g^, a) on the coupling arises from the dependence of the variational 
parameters on the coupling. The perturbative amplitude up to one- loop order is just given by 

Af[*-(g2) = const.' (l + ^\ReE{q^,0) + ilmE{q^,0)]\ . (64) 
I 27r L J J 

In the function ^(g^, a) , as opposed to the complete amplitude, there are in fact only quantita- 
tive, but not qualitative, differences between the perturbative and the variational calculations. 
For example, the imaginary parts of ImS(g2, 0) and ImS(g2, a = 0.2) both rise rapidly from 
at threshold to less than 20 just above threshold and then slowly decrease. At a = 0.8 the same 
behavior takes place, except this time the maximum of the variational calculation is a little bit 
more than twice as large. The origin of the qualitative differences between Eqs. (|63D and (|6^ 
are due to the fact that even though the coefficient a/27r is not particularly large (for exam- 
ple, it is only 0.11 for a = 0.8) - so one would expect leading order perturbation theory to be 
quite good - the numerical value of the integral ^(g^, a) more than compensates for it and 
makes the perturbative result meaningless. In particular, it is clear from Eq. ( |63D that the real 
(imaginary) parts of the amplitude will have zeros when ^ImS(g2,a) is 7r/2, 37r/2, ... (0, vr, 
...) while the perturbative result in Eq. (|64[) will not. Keeping higher terms in the perturbative 
expansion may in fact also give rise to zeros, but it is clear that, with ^ImS(g^, a) numerically 
approaching 5 at the peak for a = 0.8, to obtain any reasonable numerical agreement between 
the perturbative and the variational results one would have to expand the perturbation theory 
to a very high order indeed. 

6 Conclusion and Outlook 

The main result of this paper is a compact expression for the truncated Green function of an 
on-shell scalar nucleon interacting with an arbitrary number of mesons (Eq. (1^)), calculated 
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Figure 7: The imaginary part of the amphtude as a function of the momentum transfer squared, 
for from top to bottom, a = 0.2 (a), 0.4 (b) and 0.8 (c) The sohd hne corresponds to the use 
of the 'improved' trial action, while the dotted hne is the perturbative result 
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at lowest non-trivial order in the variational approach of Ref. |]14|. A welcome feature of this 
non-perturbative approach to the relativistic field theoretic strong coupling problem is that 
the renormalization and the truncation of this general Green function may be done analyt- 
ically and in a self consistent manner. Furthermore, by construction the expansion of the 
Green function in the coupling constant agrees exactly with the result of renormalized 1-loop 
perturbation theory. Clearly, this assures that the small coupling limit of the calculation is 
correct numerically. This was an important characteristic in the original application of this 
variational technique - numerical agreement with both the weak and strong coupling limits was 
the reason why the famous application to the polaron problem was so successful. Furthermore, 
in the present case it has the additional consequence that the analytical structure of variational 
Green functions is at least as rich as that of 1-loop perturbation theory: cuts due to particle 
production thresholds opening up as well as poles due to particle exchange in the intermediate 
state are both contained in 6*2,71- Also, of course, the technique is fully covariant; the Green 
functions may be used for any process, be it scattering or particle production, and indeed are 
even valid outside the physical regions. 

Dependence on the coupling constant of the theory enters the Green function both explicitly, 
as well as implicitly through the variational parameters. The explicit dependence allows one to 
identify contributions to individual Feynman diagrams. Indeed, one finds that the expression 
for the Green function sums up, in an approximate way, Feynman diagrams up to infinite order 
in perturbation theory (the untruncated Green function in fact contains parts of all possible 
Feynman diagrams). 

Most of the features mentioned above are independent of the functional form of the trial 
action. The analytic structure of the Green function G2,n; however, does seem to depend on this 
functional form. For the present we have been satisfied with a particular (but still very general) 
choice for the trial action (Eqs. (PB| ) and (^^)) which does allow the non-trivial truncation of 
the external dressed nucleon propagators to be made. Indeed, it results in the same variational 
equations which have already been found in Refs. |TB[ for the propagator, independent 
of the number of external mesons. This has the great practical advantage in that it means 
that the determination of the variational parameters only has to be done once. It would be 
interesting to see if this is a general feature or whether it is specific to the class of trial actions 
considered here. 

In the second part of this paper we investigated the case of meson absorption or emission and 
nucleon pair-production, i.e. n = 1, numerically. The much richer case of particle scattering, 
= 2, is considered in a following paper. After a trivial rotation to Minkowski space, the 
general expression for G2,n for n = 1 immediately yields the amplitude in the physical region 
for meson absorption/emission (i.e. < 0) as well as in the unphysical region < < 4M^. 
The representation is not defined for > 4M^, but may be analytically continued into this 
region. One finds a cut starting at the pair production threshold. At this point it should 
be remembered that we are working in the quenched approximation for the nucleons, so the 
numerical results in this region will differ to the unqenched ones. Comparison to other quenched 
calculations, such as perturbation theory, do of course remain meaningful. 

We have numerically investigated the vertex function as a function of the momentum trans- 
fer and the strength of the coupling. A strong coupling limit of the Wick-Cutkosky model does 
not exist because of its instability, so a numerical appraisal of the variational method for 
very large a has to await its application to a different field theory. Nevertheless certain non- 
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perturbative properties can already be addressed in the range of coupling constants available in 
the variational treatment of the Wick-Cutkosky model. Among these are the critical coupling 
and the width of scalar nucleon for even larger coupling constants which have been extracted 
from the propagator in Ref. |jl5|. The vertex function evaluated in the present paper provides 
another example of large nonperturbative effects: even for the relatively 'small' couplings al- 
lowed in the Wick-Cutkosky model the results differ markedly from what one would obtain in 
perturbation theory, both quantitatively and qualitatively. The origin of this disagreement may 
be traced back to the fact that although the coupling is indeed relatively small, the integrals 
multiplying the coupling are in fact very large, in particular near threshold - both in pertur- 
bation theory as well as in the variational calculation. The lesson to be learned, well known 
of course, is that smallness of the coupling is not a good criterion to determine whether or not 
to trust a perturbative result. In the present case, (1-loop) perturbation theory is inadequate 
even for values of the dimensionless coupling constant as small as 0.4 . 

In conclusion, we think that the variational method described here is a rather promising 
technique for systematically calculating Green functions in problems where perturbation theory 
is inadequate. We have shown that it is possible to truncate Green functions with one nucleon 
and an arbitrary number of external mesons in a manner which is consistent with the variational 
ansatz. It remains an interesting problem to extend this treatment to several external nucleons 
and to study bound-state problems. In addition, it seems worthwhile to find out how this 
approach will fare in realistic field theories. Attempts to apply it in QED are currently under 
way. 
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Appendix : Perturbative Limit 



It is a straightforward, but nevertheless useful, exercise to verify that the truncated (2 + n)- 
point function in Eq. (^Of ) has the correct perturbative limit. We shall sketch the derivation of 
the perturbative result in this appendix. It suffices to consider a single perturbative diagram 
^2'n(P)P'i of the form shown in Fig. |]. The complete amplitude is then obtained by 
a summation over all a and b as well as a summation over all permutations of the external 
meson momenta gj. To start off with, let us consider the case where a > b + 1, i.e. a diagram 
which has at least one external meson coupling to the nucleon while the internal meson is being 
exchanged. For this type of diagram one may replace, at this order, the unrenormalized nucleon 
mass Mq by the physical mass M. The perturbative truncated Green function corresponding 
to this situation is given by 



G2i{P,P;{(li}) = ir 



1 

2^ 



29 



b+l 



2^? 



.'i=a 



(27r)4 P + m2 



n '-^ 



(A.i; 



It is customary at this stage to rewrite the denominator in terms of integrals involving Feynman 
parameters. For the purposes of comparison to the variational calculation this is, however, not 
very useful as the variational calculation gives rise to an expression in terms of integrals over 
proper time. Although in general there will be a relation between these proper times and the 
Feynman parameters, it will be non-trivial. 

Rather, it is easier to proceed by exponentiating all denominators in the usual way: 



1 



1 



p2 + M2 2 



(A.2) 



and then performing the gaussian integration over the internal meson momentum k. If one 
exponentiates all propagators in this way one obtains a product of n + 2 integrals, the same 
number which one obtains if one expands the exponential appearing in the truncated Green's 
function in Eq. (|49|). Indeed, the n + 2 integration parameters of these two expressions may be 
identified individually if one chooses them in the following way: A nucleon propagator involving 
pf or {pi — kY between external mesons should just be exponentiated with the parameter a^. 
The two pairs of nucleon propagators between an external meson and the internal meson need 
to be exponentiated a little bit differently - one needs to recombine the integration parameters 
in the following way: 



dx dy exp 
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-^ipl + M') 



exp 



-^iiPa - kf + M') 
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exp 
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{{Pa-kf-pD 



(A.3) 
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and similarly for the two propagators involving ph+i- Here we have used the transformations 
aa = X + y and ti = Tq — x, the time being defined to be the same sum of a[s as used 
in the main text. The reason for this particular transformation should be clear from Fig. 0. 
Finally, the integration parameter associated with the internal meson should be transformed 
in the same way as the parameter u in the main text. 
Using these definitions, one readily obtains 



.i=2 



Oil 



-^{m'+pI) 

"1 du 



\ dti dt2 — e my/a, -= — , u] . (A. 4) 
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To complete the perturbative calculation of the {n + 2)-point function, we still need to 
consider the trivial case where a = b + 1, i.e. the class of diagrams where the internal meson 
loop just dresses an internal nucleon propagator. This is the only type of diagram which 
diverges and hence needs to be renormalized. For this class of diagrams one obtains, together 
with the tree level term and using the same type of exponentiation as in the derivation of 
Eq. (KM, that 
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(A.5) 



Finally then, we obtain the complete perturbative result for the truncated (2+n)-point function 
by summing over all permutations of the external momenta as well as over a and h: 

n n—1 n+1 

GIT{P,P';{<1}) = T.Gi%,p' ;{<!}) + E E G^^f{p,p';{q}) 
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Remembering that AiE) and A have perturbative expansions of the form 1 + 0{g^) (and 
hence /i^(a) = a + 0{g^)), it is easily checked that this agrees with the corresponding result 
of the variational calculation given in Eq. (^91) . 
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